WAVES, CORIOLIS FORCE AND THE DYNAMO EFFECT 
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Dynamo activity caused by waves in a rotating magneto-plasma is investigated. In astrophysical 
environments such as accretion disks and at sufficiently small spatial scales, the Hall effect is likely to 
play an important role. It is shown that a combination of the Coriolis force and Hall effect can produce a 
finite a-effect by generating net helicity in the small scales. The shear/ion-cyclotron normal mode of the 
Hall plasma is the dominant contributor to the dynamo action for short scale motions. 
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ABSTRACT 



1. INTRODUCTION 



stably stratified atmosphere (Spruit 2002), are some of 
the examples reported in the literature. 



In astrophysical objects, large scale magnetic fields 
are thought to be generated by helical turbulence (the 
so-called a-effect) and differential rotation (the il- 
effect) (see Meneguzzi, Frisch, & Pouquet (1981); 
Brandenburg (2001)). However, note that a large scale 
magnetic field can in some cases be generated with- 
out fully helical turbulence or a net a-effect. Within 
the mean field approximation, there are physical ef- 
fects which contribute to the mean electromotive force 
even when the a coefficient is zero. A shear turbulent 



Helicity is naturally imparted to a rotating fluid 
by the Coriolis force (Moffatt 1970a,b, 1972, 1978). 
However, at sufficiently small scales, where the 
Rossby number 



is larger than unity (Uq and L are characteristic veloc- 
ities and lengths, and is the rotation rate), the Cori- 
olis force and the resulting induced helicity might be- 
come negligibly small. Nonetheless, it is worth noting 
that in helical turbulent flows, the kinetic helicity de- 
velops a direct cascade along with the energy (Cheng 
et al. 2003; Gomez & Mininni 2004). Therefore, some 
small scale turbulent flows can still be helical, even 
though the source of helicity remains at much larger 
spatial scales. 



Rs 



2L n 



u 



(1) 



flow (Urpin 2002), the lu x j term (Riidler 1969; Gep- 
pert & Rheinhardt 2002), or magnetic instabilities in a 
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The Hall effect introduces a definite handedness or 
helicity on small scale fluid motions (Wardle (1999), 
Balbus & Terquem (2001)), since the mirror symme- 
try in the induction equation is broken (see Eqn (3) 
below). Therefore, one should expect a net a-effect in 
a Hall plasma. The Hall effect becomes relevant when- 
ever the Hall length scale 



U 



(2) 



is larger than the dissipation scale, a category to 
which several objects of astrophysical interest belong 
(Mininni, Gomez, & Mahajan 2002a, 2003b). Here, 
Ua is the characteristic Alfvenic speed, c is the speed 
of light, and uj pi is the ion plasma frequency. 

In §2 we write down the Hall MHD equations. The 
normal modes sustained in this system are derived and 
listed in §3. In §4 we briefly summarize the role of 
normal mode fluctuations on MHD dynamos, and in 
§5 this concept is extended to Hall MHD. In §6 we 
show the effect of rotation on the a-effect. The main 
results of the present work are summarized in §7. 

2. THE HALL-MHD SYSTEM 

The dynamics of ideal and incompressible fully 
ionized plasmas in a rotating frame, is described by 
the induction equation (modified by the addition of the 
Hall current) and the Navier-Stokes equation, 



dB 

dt 
dU 

~dt 



= Vx[([/-AVxB)xB] (3) 
= -([/• V)U-2nxU + (X7xB) xB 



- v(P- |fJxrf) 
with the constraints 



V U = V • B = . 



(4) 



(5) 



These equations are known as the Hall-MHD equa- 
tions. The magnetic field is expressed in velocity units, 
i.e. B — B[Gauss](4Trp)^ 1 ^ 2 , where p is the con- 
stant mass density. The quantity P is the gas pressure 
divided by the constant mass density. We define a di- 
mensionless number e to measure the relative strength 
of the Hall effect, 



e = 



To' 



(6) 



where Lq is the characteristic length scale of the sys- 
tem. If we choose Uq = Ua as the characteristic veloc- 
ity, A reduces to the ion skin depth. Note that equation 
(2) is valid for a fully ionized plasma. 

3. WAVES IN HALL-MHD 

We will study waves in the ideal and incompressible 
Hall-MHD system in a rotating frame. To linearize 
the Hall-MHD equations around a static and uniform 
magnetic field Bq, we write 



B = B Q + b 

U = u, 



(7) 
(8) 



where u,b ~ exp(ifc • r — iut). These substitutions 
convert (3) and (4) into the closed set : 

ub = kx[(u-i\kxb) xB ] (9) 
-urn = 2iftxu + (kxb) xB - kPtot (10) 

where the total effective pressure P tot — P — | f2 X r | 2 . 

The elimination of P to t in incompressible flows is 
arranged by projecting onto the plane perpendicular to 
fe . One finally obtains 

-ub = (fe • B a ) (u - iXkxb) (11) 
-lou = V[2inxu+ (kxb) XB ] , (12) 

where Pjj = 5ij — kikjk~ 2 is the projector operator. 

Without loss of generality we choose k in the z- 
direction, i.e. k = kz, while ft and B can be oriented 
arbitrarily. Equations (11) and (12), then, reduce to 

(cjI + i\k 2 B z A) b = -kB z lu (13) 

(loI — 2iQ z A) u = —kB z Ib . (14) 

Here, 



I : 



and 




(15) 



(16) 



Note that only the z-components of ft and B are 
relevant (i.e. the components in the direction of the 
vector fc). The orthonormal base for the antisymmetric 
operator A is given by 



l±> 




(17) 
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and satisfies A |±) = ±i |±). 

The dispersion relation can be obtained from Equa- 
tions (13) and (14), 

lu 2 - k 2 (B 2 + 2XQ Z B Z ) + au (2fl z - \k 2 B z ) = , 

(18) 

where a = ±1 and the eigenvectors can be written as 
\a). 

This dispersion relationship is quite general and 
includes several well known waves of a magnetized 
plasma in the corresponding asymptotic limits. When 
both Q, z = (no rotation) and A = (negligible Hall 
current, i.e. the MHD limit) we obtain uj 2 — k 2 B 2 = 
which corresponds to Alfven waves. In the incom- 
pressible limit, the shear and the compressional waves 
are degenerate and indistinguishable. When A = 0, 
equation (18) reduces to lu 2 + 2auuVL z — k 2 B 2 = 0, 
and we obtain the inertial waves first described in the 
absence of an external magnetic field B z by Moffatt 
(1970a,b, 1972). Finally, when fi z =0we obtain two 
branches of circularly polarized waves. The right hand 
polarized branch, corresponds to whistlers, with a fre- 
quency growing like uj ~ k 2 at large wavenumbers. 
These waves are the high k limit of the compressional 
branch. The left handed branch is the standard shear 
wave whose frequency approaches the ion-cyclotron 
frequency asymptotically. 

Figure 1 shows these asymptotic cases in detail. 
The general dispersion relationship given in equation 
(18), can be cast in dimensionless units using Lq and 
B z = Ua as typical longitude and velocity, 



k 2 (l + pe) + auu(p - ek 2 ) = , 



where 



P 



2n z L 



B, 



(19) 



(20) 



can be interpreted as the inverse of the Rossby number 
given in Eq. (1). The MHD limit in Figure 1 cor- 
responds to a small neighborhood around the origin, 
where there is only a transition between inertial waves 
and Alfven waves. 

Figure 2 shows the phase speed as a function 
of wavenumber for the two positive branches given 
by equation (19). There are also two other identi- 
cal branches with negative frequencies (not shown). 
These waves behave approximately as Alfven waves 
only in the wavenumber regions where they are non- 
dispersive, i.e. when the curves become horizontal. 



4. MHD WAVES AND THE DYNAMO EFFECT 

Before embarking on the Hall-MHD case, it is in- 
structive to recall some previous results relating to the 
dynamo process induced by MHD waves ( readily ob- 
tained from our expressions). In the limit f2 — > 
and A — > 0, Alfven waves are known to quench the 
a-effect (Gruzinov & Diamond 1994). The MHD a- 
coefficient (Pouquet, Frisch, & Leorat 1976; Black- 
man & Field 1999) is given by 

a =I(-( M .VxM} + (5-Vx6}) , (21) 

where u and b are respectively small scale velocity 
and magnetic fields, unaffected by the presence of a 
large scale magnetic field (see details in Blackman & 
Field (1999)), and the coefficient r is a typical correla- 
tion time for the turbulent small scale motions. A pure 
Alfvenic state satisfies u = ±b and therefore a = 0. 
This is to be expected, since all nonlinear terms cancel 
exactly for Alvenic states, and therefore no transporta- 
tion coefficients can arise in that case. On the other 
hand, when f2 ^ the Coriolis force is expected to 
inject helicity to the fluid, and therefore a net a-effect 
arises (Moffatt 1978). A detailed discussion can also 
be found in Moffatt (1970a,b, 1972) 

5. HALL-MHD NORMAL MODES AND THE 
a-EFFECT 

The first studies on the impact of Hall currents on 
dynamo action (Helmis 1968; Helmis 1971) were car- 
ried out using mean field theory and the first-order 
smoothing approximation (Krause & Radler 1980). 
Helmis obtained decreasing dynamo action as the 
strength of the Hall terms increased. Recently, the 
relevance of the Hall effect on dynamo activity was 
confirmed experimentally by Ding et al. (2004). 

In Mininni, Gomez, & Mahajan (2002a) it was 
shown that the expression for the a-effect in the pres- 
ence of the Hall effect is modified according to 



a = -(-(tt e -Vxtt e ) + (b-Vx6)- 



-A (6- VxVx« e )) 



(22) 



where u e = u — AVxb is the small scale electron 
flow velocity. This general expression differs from the 
MHD result (equation (21)) in two ways: it replaces 
the kinetic helicity by the helicity of the electron flow, 
and it contains an extra term due to the Hall current 
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in the microscale. A nontrivial consequence of the lat- 
ter is that, while the original a coefficient of Pouquet, 
Frisch, & Leorat (1976) is zero for a pure Alfvenic 
state u = ab (Gruzinov & Diamond 1994), the one 
corresponding to equation (22) is not. In Mininni, 
Gomez, & Mahajan (2003b) the impact of the Hall ef- 
fect in helical turbulent dynamos was studied in direct 
numerical simulations. As mentioned in the previous 
section, the required helicity is naturally introduced in 
a fluid in a rotating body by the Coriolis force. 

Two questions arise. The first one is better posed 
and answered for a non-rotating plasma. We just 
showed that the a derived in MHD goes to zero for 
the pure Alfvenic state , which is an eigen-state of 
MHD. We also claimed that for the same pure eigen- 
state of MHD, the a derived in Hall MHD does not 
vanish. Surely that is an unwarranted mixing of differ- 
ent worlds. What we should, instead, calculate is the 
value of a derived in Hall MHD for the corresponding 
normal modes of Hall MHD. This task is performed at 
the end of this section. 

The second question concerns rotation and the 
Coriolis force which is known to inject helicity at large 
scales (see equation [1] and Figure 1). At small scales 
(large k), however, the Coriolis force is not considered 
to be relevant (motions are expected to be essentially 
non-helical) posing a serious restriction on the gener- 
ation of magnetic fields in astrophysics. Can we find a 
source of helicity acting at small scales, and will this 
source have much to do with rotation? 

We show that the Hall effect naturally introduces 
helicity at small scales, which is precisely the region 
where this effect is stronger (in agreement with Wardle 
(1999) and also Balbus & Terquem (2001)). However, 
no net a-effect is generated by these microscale mo- 
tions, unless there is also a net rotation of the system; 
a combination of rotation and Hall effect is needed for 
dynamo action. 

As is shown in equation (17), the general solu- 
tions of the linearized Hall-MHD equations in a ro- 
tating frame are right-handed or left-handed polarized 
waves. To investigate the effect of Hall currents at 
small scales, we first concentrate on a non-rotating 
plasma, Q z = 0. For large k, the dispersion relation- 
ship (18) reduces to 

u 2 - k 2 B 2 z - a\uk 2 B z = , (23) 

which coupled to Eq.(l 1), yields 

b=-^u e . (24) 



Inserting Eq. (24) into equation (22), and invoking that 
the electron vorticity is Vxit e = aku e , we obtain 

r fu 2 - k 2 B 2 - a\ujk 2 B z \ e . ^ 

a = -- f z — 2 J (u e • Vxit e ) (25) 

which, according to the dispersion relation Eq. (23), 
corresponds to a = 0. This result of a zero a-effect for 
Hall MHD normal modes is an expected and natural 
extension of the similar result obtained for Alfvenic 
states in MHD. However, we must bear in mind that 
these results are derived for non-rotating systems. 

6. CORIOLIS FORCE AND THE a-EFFECT 

In rotating objects, the closure calculation leading 
to either equation (22) (in the MHD limit) or equation 
(21) needs to be revised, since the Coriolis force (see 
equation (4)) was not included in those calculations. 
The starting point is the so-called "reduced smoothing 
approximation" (RS A) proposed by Blackman & Field 
(1999) (see also Mininni, Gomez, & Mahajan (2003a) 
for a derivation which includes the Hall effect). Fol- 
lowing RSA, we decompose the magnetic and velocity 
fields as 

B = B + b + b° (26) 

U = U + u + u° (27) 

where the overbar denotes spatially or statistically av- 
eraged large-scale perturbations. The small scale fields 
b°, u° are solutions of Eqs. (3)-(4) in the absence 
of large scale fields, and b, u are anisotropic correc- 
tions to the small scale fields, caused by the presence 
of B, U. The net effect of small scale fluctuations on 
the large scale dynamics, is given by an electromotive 
force 

£ = {u°' e x b + u e x b°) (28) 

acting on equation (4), where u e = u — AV x b and 

u o,e = u o _x\7 x b°. 

From equations (3)-(4) for the evolution of small 
scale fields (assuming that U w in an appropriate 
frame of reference, and under the reduced smoothing 
approximation (Blackman & Field 1999)), we obtain 

d t b ~ (B ■ V) «°< e , (29) 

[dt + 2I-nx]u~ (B- V)6° . (30) 
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The time derivative indicated in equations (29)-(30) 
is usually approximated by dt ~ 1/r, where r is a cor- 
relation time for the microscale motions. The operator 
on the left-hand side of equation (30) becomes, 



T -i 



+ 21 -fix 



(31) 



whose inverse is 



T = 



i + (2n Z T) 



(I + 2^tA) , (32) 



Replacing equations (29)-(30) on equation (28) and 
using equation (32), we obtain 



£ = 



TokB, 



-u 2 + 



k 2 B\ 



o 2 V ~ ' l + (2ft z r) 2 
+ Xak 2 B x u) (|w°' e | 2 ) . 



+ 



(33) 



The expression for the electromotive force given in 
equation (33) corresponds to an anisotropic tensor 



a« = -^f <|u°' e | 2 )[2fi z (Afc 2 I? 

2 

.2 d2 



+k z B 2 



(2» z r) 2 
1 + (20 z r) 5 



k^kj 



(34) 



which for small scale fluctuations given by a pure 
mode of wavenumber fc, produce an electromotive 
force £ parallel to k, regardless of the orientation of B 
and fi. It seems reasonable to assume that the corre- 
lation time t is much smaller than the rotation period. 
Therefore, in the asymptotic limit fir <C 1, 



2fl z rk 



Oli 



<|«°-«| a > 



\ak 2 B 7 



- 1 



k% kj 



(35) 

Note that when a rotation field is present, ay is in gen- 
eral non-zero. The rotation field therefore provides 
a source of energy that is responsible for the net a- 
effect. Equation (35) confirms that for non-rotating 
objects, the net alpha effect generated by a background 
of small scale normal modes is exactly zero. Notwith- 
standing, sources of kinetic helicity other than rota- 
tion have been considered (Mininni, Gomez, & Maha- 
jan 2003a,b), to assess their efficiency in driving large 
scale dynamos. 

For large wavenumbers, the dispersion relationship 
given by equation (18) has two limits: one is the so 
called compressional/whistler branch ui w ±Xk 2 B z , 
while the other is the shear/cyclotron branch with 



ui s=s ± B z /\. For the former, the a-effect is asymp- 
totically small. For the shear/cyclotron branch, how- 
ever, otij oc fc 3 at large wavenumbers implying strong 
dynamo action. Therefore, once the rotation breaks 
the mirror symmetry, the shear/ion-cyclotron modes ( 
which failed to produce dynamo action without rota- 
tion) are, indeed, able to provide a net a-effect at mi- 
croscopic scales. We believe that this is a very impor- 
tant result for turbulent dynamo theories. In Figure 3 
we show the trace of the tensor a, 3 as a function of fc, 
for different values of p and e. Note that in the MHD 
case an drops to zero at small scales (large fc), while 
in the Hall-MHD case it does not. 

7. DISCUSSION 

We have shown that the Hall effect in conjunction 
with fluid rotation can generate helicity at small scales 
(i.e. produce small scale helical motions) leading to 
a net a-effect through the agency of the shear/ion- 
cyclotron normal mode of the plasma. This finding 
can be of considerable importance to the existence of 
large scale dynamo action in a variety of astrophysi- 
cal objects. Our results are quite consistent with pre- 
vious results obtained in the study of instabilities in 
accretion disks. Wardle (1999), Balbus & Terquem 
(2001), and Sano & Stone (2002) showed that the 
magneto-rotational instability can be either enhanced 
or quenched by the Hall effect depending on the ori- 
entation of 1~2 and Bo, which is just a manifestation 
of the handedness introduced by the Hall effect. In a 
future work the detailed investigation of this mode of 
dynamo action will be carried out through direct nu- 
merical simulations. 
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Fig. 2. — Phase speed vs. wavenumber for the normal 
modes arising from the normalized dispersion relation- 
ship given in equation (19). Each frame corresponds to 
a different value of e, while different traces correspond 
top = 0.01 (thick) and p = 0.1, 1, 10, 100. 
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Fig. 3. — The trace of aij vs. k in arbitrary units (see 
equation (35) for the two positive branches of the dis- 
persion relationship. Each frame corresponds to a dif- 
ferent value of the dimensionless rotation speed p (in- 
dicated). The dotted curves correspond to e = 0, the 
thin full curves to e = 0.1, and the thick full curves to 
e = 1. 



